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Abstract. We present sub- millimetre images of the spi- 
ral galaxy NGC 6946 at 450 and 850 /mi. The observa- 
tions have been carried out using the scan/map imag- 
ing mode of the Sub-millimetre Common User Bolometer 
Array (SCUBA) at the James Clark Maxwell Telescope 
(JCMT). The comparison with a 12 CO (2-1) line emission 
map from IRAM and with an U-band image, clearly shows 
that the 850 /im emission is associated with molecular gas 
and bright star-forming regions, tracing the spiral arms. 
We place an upper limit of 0.7 for the V-band optical 
depth of dust associated with a diffuse atomic component. 
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1. Introduction 

Until recently,the main source of Far-Infrared (FIR) data 
for spiral galaxies has come from the IRAS satellite. The 
availability of instruments capable of observing at A > 
100 /im, like those onboard the satellite ISO, has allowed 
the detection of large amounts of cold dust, much colder 
than IRAS was able to detect. Spiral galaxies are found 
to have a dust content similar to the one in the Galaxy 



(Alton et al. |l998aj ). 

Cold temperatures (T<20K) can be reached by dif- 
fuse dust heated by the general interstellar radiation field, 
while du st clo se to star-forming region is hotter (T>50K; 
Whittet 1992 ). Since diffuse dust is the main contributor 
to the internal extinction in a galaxy, observations of cold 
dust help to trace its opacity. Additionally, ISO obser- 
vations have suggested that the dust distribution is more 
extended than the stellar disk (Alton et al. 1998a; Bianchi, 
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Davies & Alton 1999b). If this is confirmed, observation 
of the distant universe may be severely biased, because of 
the large cross section of the dust disks. 

Unfortunately, the poor resolution (~ 2') of the current 
FIR images does not allow detailed studies of the spatial 
distribution of dust. Higher resolution can be achieved in 
the sub-mm, but a high sensitivity is required because of 
the fainter dust emission. High sensitivity and resolution 
are both characteristics of the recently developed SCUBA 
sub-mm camera. Only a few large nearby galaxies have 
been observed with SCUBA, notably the highly inclined 
galaxy NGC 7331 (Bianchi et al. 
galaxy NGC 891 (Alton et al. 



1998 ) and the e dge-o n 
1998b| ; Israel et al. |l999|) . 



The observed dust emission is found to correlate well with 
the molecular gas phase, dominant in the centre. How- 
ever, a dust component associated with the atomic gas is 
needed to explain the dust and gas column density at large 
galactocentri c dist ance along the major axis of NCG 891 
(Alton et al. |1999| ). 

In this Letter we present SCUBA observations of the 
face-on galaxy NGC 6946. Because the galaxy is larger 
than the camera field of view, images have been produced 
with the scan-mapping technique, chopping within the 
observed field. The observation and the data reduction 
needed to restore the source signal are described in the 
next section. The description and the discussion of the 
results are given in Section 3. 

2. Observations and data reduction 

NGC 6946 was observed at 450 /im and 850 /im, during 
April 10, 11 and June 17, 18, 19, 20 1998. 

SCUBA consists of two bolometer arrays of 91 ele- 
ments optimised to observe at 450 /im and 37 elements 
optimised at 850 /im, covering a field of view of about 2.3 



arcmin (Holland et al. 1999). The camera, mounted on the 
Nasmyth focus of the telescope can be used simultaneously 
at both wavelengths, by means of a dichroic beamsplitter. 

In the scan-map mode, the telescope scans the source 
at a rate of 24 arcsec per second, along specific angles 
to ensure a fully sampled map. Meanwhile the secondary 
chops with a frequency of 7.8 Hz within the observed field. 
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While this ensures a correct subtraction of the sky back- 
ground, the resulting maps unfortunately have the profile 
of the source convolved with the chop. The profile of the 
source is restored deconvolving the chop from the observed 
map by mean of Fourier Transform (FT) analysis. 

Scan-maps of NGC 6946 presented here are fully sam- 
pled over an area of 8'x8'. Each set of observations con- 
sisted of six scans, with different chop configurations: chop 
throws of 20", 30" and 65" along RA and Dec are needed 
to retrieve the final image. Data have been reduced us- 
ing th e STARLINK package SURF (Jenness & Lightfoot 
1999). Images were first flat-fielded to correct for differ- 



ent sensitivities of the bolometers. Noisy bolometers were 
masked and spikes from transient detections removed by 
applying a 5-tr clip. A correction for atmospheric extinc- 
tion was applied, using measures of the atmosphere opac- 
ities taken several times during the nights of observation. 
Zenith optical depth varied during the six nights, with 
7450 = 0.4 — 2.5 and Tgso = 0.1 — 0.5. The 450 //m opacity 
on the last night was too high (r > 3) for the source to 
be detected and therefore the relative maps were not used 
for this wavelength. Because of the chopping in the source 
field, each bolometer sees a different background: a base- 
line, estimated from a linear interpolation at the edges of 
the scan, has been subtracted from each bolometer. 

Sky fluctuations were derived from the time sequence 
of observations for each bolometer, after the subtraction of 
a model of the source, obtained from the data themselves. 
The images have then been corrected by subtacting the 
systematic sky variations from each bolometer. 

Data taken with the same chop configuration were re- 
binned together into a map in an equatorial coordinate 
frame, to increase the signal to noise. Six maps with 3" 
pixels were finally obtained for each wavelength, combin- 
ing 33 and 25 observations, at 850 and 450 /im, respec- 
tively. In each of the six maps the signal from the source is 
convolved with a different chop function. The final decon- 
volved image is retrieved using the Emerson II technique 
(Holland et al. |1999[ Jenness, Lightfoot & Holland |l99§) . 
Essentially, for each rebinned image, the FT of the source 
is derived by dividing the FT of the map by the FT of 
the chop function, a simple sine-wave. Since the division 
boosts up the noise near the zeros of the sine-wave, dif- 
ferent chop configuration are used. For the chosen chop 
throws, the FT of the chop functions do not have coinci- 
dent zeros, apart from the zero frequency. A smoother FT 
of source can therefore be obtained, and the final image is 
retrieved by the applying an inverse FT. 

Unfortunately the deconvolution introduces artifacts 
in the images, like a curved sky background. This may be 
due to residual, uncorrected, sky fluctuation at frequencies 
close to zero, where all the chops FT goes to zero. Work 
to solve this problem is ongoing (Jenness, private commu- 
nication). To enhance the contrast between the sky and 
the source, we have modelled a curved surface from the 
images, masking all the regions were the signal was evi- 



dently coming from the galaxy. The surface has been then 
subtracted from the image. 

Calibration was achieved from scan-maps of Uranus, 
that were reduced in the same way as the galaxy. Inte- 
grated flux densities of Uranus were derived, for each ob- 
serving period, using the STA RLIN K package FLUXES 
(Privett, Jenness & Matthews |1998| ) for JCMT planetary 
fluxes. Comparing data for each night we derived a rela- 
tive error in calibration of 8 per cent and 17 per cent, for 
850 /im and 450 /im respectively. From the planet profile, 
the beam size was estimated: FWHMs of 15.2" and 8.7" 
were measured for the beam at 850 and 450 /mi, respec- 
tively. To increase the signal to noise, the 850 /im image 
has been smoothed with a gaussian of 9" thus degrading 
the beam to a FWHM of 17.7" The 450 /im image has 
been smoothed to the same resolution as for the 850 /im 
one, to facilitate the comparison between features present 
in both. The sky a in the smoothed images is 3.3 mJy 
beam -1 at 850 /im and 22 mJy beam" 1 at 450 /im. 

The final images, after removing the curved back- 
ground and smoothing are presented in Fig. [IJ. For each 
wavelength, the grey scale shows all the features >l-er, 
while contours starts at 3-cr and have steps of 3-cr. 

3. Results and discussion 

The 850 /xm image shows a bright nucleus and several 
features that clearly trace the spiral arms (in Fig. [I] the 
sub-mm contours are ov erlayed to a U-band image of the 



galaxy (Trewhella 1998)) As already seen in optical im 



ages (Tacconi & Young |l990J) , the spiral arms originating 
in the northeast quadrant are more pronounced than the 
others, where only regions with bright HII regions have 
detectable emission in the sub-mm. The 850 /im image 
presents a striking similarity to the 12 CO(2-l) emission 
map in Sauty, Gerin & Casoli ( |l998j ), observed with the 
IRAM 30m radiotelescope with a comparable resolution 
(13"). The image of molecular line emission is also shown 
in Fig. [|, with the 850 /im contour overlayed. The similar- 
ity with the sub-mm image is hardly surprising, since the 
molecular gas is the dominant component of the ISM over 



the optical disk of NGC 6946 (Tacconi & Young |1986J). 
The nucleus is elongated in the direction north-south, as 
observed for the central bar of molecular gas (Ishizuki et 
al. [1990|; Regan & Vogcl |1995[). 



Emission associated with a more diffuse atomic gas 
component cannot be detected, for several reasons. First 
of all the face-on inclination of the galaxy: since dust is 
optically thin to its own emission a faint component can 
be observed only if the dust column density is large. This 
is the case for the high inclination galaxies NGC 7331 



(Bianchi et al. 1998) and NGC 891 (Alton et al. 1998b) 



where higher signal to noise were obtained coadding a 
smaller number of observations. The large face-on galaxy 
M51 has been observed using the scan-map mode and con- 
firms the necessity of long integrations (Tilanus, private 
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Fig. 1. Sub-mm images of NGC 6946, at 850 /im (top-left) and 450 /im (top-right). Grey scales show features 1-a 
above the sky, while contours starts at 3-ct and have steps of 3-cr. Both images have a beam size FWHM=17.7" An 
area of 10'xlO' is displayed, but only the central 8'x8' are fully sampled. North is on top, Eas t on t he left. A U-band 
image of NGC 6946 (Trewhella 1998 ) and a 12 CO(2-l) emission map (Sauty, Gerin & Casoli 1998 ), are presented in 
the bottom-left and bottom-right panels, respectively. 850 /im contours are overlayed to both optical and line e missio n 
images. The centre of the galaxy is in RA=20' 1 3 4 m 5 2!3 and Dec=60° 09' 14'.'21 (J2000; De Vaucouleurs et al. |l99l|) . 
The scale is 27 pc arcsec" 1 (D= 5.5 Mpc; Tully |l988[). 



communication). Furthermore, chopping inside the source 
field removes not only the emission from the sky but also 
from possible components with a shallow gradient: this 
may be the case for dust associated wit h the fl at HI dis- 
tribution in NGC 6946 (Tacconi & Young |l986| ). Finally, a 
faint diffuse emission could have been masked by the men- 



tioned artifacts and subtracted together with the curved 
background. 



The 450 /im image is much noisier than the 850 fim 
one, because of the larger sky emission at this wavelength. 
Only a central region of 0'75 x can be clearly detected, 
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although most of the features at a 3-a level correspond to 
regions emitting in the long wavelength image. 

The temperature from the two sub-mm fluxes can 
be measured only for the central region with significant 
450 /mi flux. Sub-mm fluxes are 1.2 Jy at 850 /im and 
9.3 Jy at 450 /urn. We checked for the contribution of the 
strong 12 CO(3-2) line emission at 346 GHz to the 850 /im 
flux using the observation of NGC 6946 centre in this 
line reported by Mauersberger et al. ( 1999 ) for a beam of 
21". C onver ting from the original units to Jansky (Braine 
et al. |1995| ) and averaging over t he 30 GHz bandwidth 
of 850 /mi filter (Matthews |1999p , a flux density of 80 
mJy/beam is derived. As pointed out by the referee, the 
pointing of the Mauersberger et al. observations was offset 
from the strongly concentrated central emission by nearly 
one beamwidth in the SW direction. Using the 12 CO(2-l) 
image as a template of the higher state emission, we cor- 
rected for the offset and derived the flux for the central 
region, larger than the beam. A total contribution of 0.6 
Jy is derived for the 346 GHz line (50% of the 850 /mi 
flux). However, this large contribution is due to the high 
density and gas temperature of the central region. In fact, 
for 850 /im fluxes on larger apertures, the contamination is 
much smaller: Israel et al. ( |1999[ ) derive a contribution of 
only 4% to the total sub-mm flux of NGC 891. Therefore, 
the derivation of cold du st temperature at large galactic 
radii (Alton et al. 1998b| ) are not severely biased. We did 
not correct the 450 /mi flux for the contribution of the 
12 CO(6- 5) line, that lies at the edge of the filter (Israel et 



N d [10 10 em- 2 ] 



al. |1999P . 



After the correction, the dust temperature of the cen- 
tral region is T=34±6 K, where the large quoted error 
comes from the calibration uncertainties. Here and in the 
following, dust temperature and masses are computed us- 
ing the emissivity law Q cm (A) derived by Bianchi, Davies 



& Alton (1999a) from observation of diffuse FIR emission 



and estimates of optical extinction in the Galaxy. For a 
wavelength dependence of the emissivity A _/3 , changing 
smoo thly from = 1 to (3 = 2 at 200 /im (Reach et al. 
1995| ) they obtain Q GX t(V)/Q om (100/im) = 2390 ± 190, 
where Q cx t(V) is the extin ction efficiency in the V-band 
(Qcxt(V) f«1.5; Casey 1991). Lacking information outside 



of the centre, a mean temperature for a larger aperture 
can be derived from the lower resolution IRAS and ISO 



images at 100 /im and 200 /tm (Alton et al. 1998a ). The 
total flux inside the B-band half light aperture (5' in di- 
ameter) is 240±40 Jy at 1 00 /mi a nd 280±40 Jy at 200 /mi 
(Bianchi, Davies & Alton 1999b| ). The temperature from 
IRAS and ISO fluxes is T=24±2 K. 

We derived a point-to-point correlation between the 
850 /tm flux and the 12 CO(2-l) line, resampling the sub- 
mm image to the same pixel size as the line emission map 
(10", roughly equivalent to both beam sizes) and using 
all positions with signals larger than 3-a in both obser- 
vations. A linear correlation is found (Fig. ^|). Assuming 
a mean dust grain radius a — 0.1/tm and mass density 



bL 



0.1 0.3 
bO/j,m surface brightness [Jy beam L J 



Fig. 2. Point-to-point correlation between the 12CO(2-l) 
line emission and the SCUBA flux, for signals larger than 
3-cr in both observations. Column density are derived as 
described in the paper. 



p=3 g cm -3 (Hildebrand 1983| ), the emissivity of Bianchi 
et al. (1999a)^ and T=24K, the dust column density and 
hence the mass along the line of sight can be easily com- 
puted. The molecular gas column density has been derived 
from the 12 CO(2-l) emission using a conversion factor ap- 
propriate for the 12 CO(1-0) emission in the general ISM in 
the Galaxy (X=1.8 1 20 citT 2 K" 1 km" 1 s; Maloney fL990|) 



and a line ratio I(2-l)/I(l-0)=0.4 (Casoli et al. |1990p . The 
slope of the linear correlation can then be converted into 
a gas-to-dust mass ratio of 170±20, a value very close to 



the local Galactic one (160; Sodroski et al. 1994). This 
confirms the association of dust with the local dominant 
phase of the galactic ISM. 

The dust content of NGC 6946 has been studied car- 
rying out an energy balance between the stellar emission 
in the optical and the FIR dust emission, through the 
help of radiative transfer models. If an exponential disk is 
used to model the dust distribution, a central face-on op- 
tical depth tv ~ 5 is needed to explain the FIR emission 
(Evans [1992^ Trewhella |l998| Bianchi et al. |l999b| ). The 
850 /tm image clearly show that the dust distribution is 
more complex, but still the column densities derived from 
the sub-mm flux support the idea of an optically thick 
dust distribution. Under the same assumption of the pre- 
vious paragraph, the diffuse component of the north-east 
spiral arms at a 3-a level corresponds to a V-band optical 
depth ry ~ 2.2. The quite high optical depth correspond- 
ing to the sky noise (ry « 0.7) shows how difficult is to 
obtain sub-mm images of dust emission in the outskirts of 
face-on galaxies, even for a high sensitivity instrument like 
SCUBA. Thus, possible extended dust distributions (Al- 
ton et al. 1998a) are better revealed through deep sub-mm 
imaging of edge-on galaxies, where the dust column den- 



1 Using the Bianchi et al. emissivity for Q om oc A 2 at any 
wavelength, results in dust column densities smaller by only 
15%. 
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sity is maximized. However, the high inclination makes 
the interpretation of the dust emission along the line of 
sight more complex. 
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